d-wave superconductivity on the checkerboard Hubbard model at weak and strong 

coupling 
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It has been argued that inhomogeneity generally can enhance superconductivity in the cuprate 
high-Tc materials. To check the validity of this claim, we study d-wave superconductivity on the 
checkerboard Hubbard model on a square lattice using the Cellular Dynamical Mean Field theory 
method with an exact diagonalization solver at zero temperature. The d-wave order parameter is 
computed for various inhomogeneity levels over the entire doping range of interest in both strong and 
weak coupling regimes. At a given doping, the size of the d-wave order parameter manifests itself 
directly in the height of the coherence peaks and hence is an appropriate measure of the strength 
of superconductivity. The weak coupling results reveal a suppression of the order parameter in the 
presence of inhomogeneity for small to intermediate hole dopings, while it is enhanced for large 
dopings. In contrast, for strong coupling there is a monotonic decrease in the maximum amplitude 
of the superconducting order parameter with inhomogeneity over the entire doping range of interest. 
Furthermore, at moderately high inhomogeneity, the system undergoes a first-order transition from 
the superconducting to the normal state in the underdoped regime. In the overdoped regime, the 
change in the value of the superconducting order parameter correlates with the height of the lowest 
energy peak in the spectral weight of antiferromagnetic spin fluctuations, confirming the connection 
between antiferromagnetic fluctuations and d-wave superconductivity found in earlier studies on the 
homogeneous case. Our results are benchmarked by comparisons with numerically exact results on 
the checkerboard Hubbard ladder. 



INTRODUCTION 



A variety of experimental probes have revealed the 
existence of inhomogeneity in the electronic spectra of 
several families of the cuprate high-T c superconductors, 
which manifests itself as spatial modulation of the charge 
or spin density. These manifestations of inhomogene- 
ity include the one-dimensional stripe ordered phase ob- 
served in neutron scattering experiments^ and the two- 
dimensional checkerboard pattern observed in optical 
spectroscopy measurements^ and STM studies on un- 
derdoped Bi-2213^ and CaNaCuOC^. Such observa- 
tions have motivated a number of theoretical scenarios 
involving inhomogeneity as the key ingredient in high-T c 
superconductivity^—. In addition, there exist a number 
of theoretical results on the Hubbard model reporting an 
enhancement in the strength of superconductivity in the 
presence of inhomogeneity^^— , while other studies^ find 
a suppression of superconductivity with inhomogeneity. 
Whether inhomogeneity in the cuprates is a friend or foe 
of superconductivity in cuprates is therefore still an open 
issue. 

Here we employ the Cellular Dynamical Mean-Field 
Theory (CDMFT) approach at zero temperature with an 
exact diagonalization solver to study d-wave supercon- 
ductivity on what is commonly called (somewhat abu- 
sively) the checkerboard Hubbard model. In contrast to 
previous quantum cluster based studies^ using Dynam- 
ical Cluster Approximation (DCA) at finite temperature 
that focused entirely on the underdoped regime, we con- 
sider the entire doping range of interest from half-filling 
up to the extreme overdoped regime and over a wider 
range of inhomogeneity. In addition, we study both the 



weak and strong coupling regimes. 

In Section II we introduce the checkerboard Hubbard 
model and describe the details of the method employed 
to study superconductivity. In Section III we present 
comparisons with previous accurate results^ that serve 
as benchmark for our approach. We present our results in 
Section IV. This is followed up in Section V by discus- 
sions and further comparisons with previous literature. 
Finally, a summary of main results and final conclusions 
appear in Section VI. 



II. METHOD 

Our starting point is the one band Hubbard model on 
a two-dimensional square lattice 
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in which electrons hop among a set of lattice sites, but 
pay an energy cost U whenever they doubly occupy the 
same site. Here i,j label lattice sites, the hopping ma- 
trix elements Uj vanish unless i,j are nearest neighbors, 
and di a annihilates an electron with spin a on site i. 
As a simple toy model for inhomogeneity based on this 
Hamiltonian, we consider a checkerboard modulation of 
the nearest neighbor hopping amplitude tij which varies 
between alternate bonds with values t and if along cither 
direction, (Figure [1]) with t = t' being the homogeneous 
case. 

We employ CDMFT, a cluster generalization of Dy- 
namical Mean Field Theory (DMFT)±£ that allows one to 
reliably study d-wave superconductivity. In CDMFT, the 
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FIG. 1. Checkerboard lattice with hopping amplitudes t 
(thick lines) and t' on alternate bonds. A four site plaquette 
can be chosen in three distinct ways namely: (1) all t bonds 
(plaquette ^4), (2) t along y(x) and t' along x(y) (plaquette 
B) and (3) all t' bonds (plaquette C) 



lattice problem is mapped to one involving a finite clus- 
ter coupled to a bath of non-interacting electrons^—. 
The local quantum correlations within the cluster are in- 
cluded exactly while longer-range correlations are treated 
using a mean-field approximation by writing down an ef- 
fective action 
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where Go is a dynamical (time dependent) Weiss field 
that describes the coupling of the cluster to the bath. 
The cluster is a four-site (2 x 2) plaquette, which has 
been used extensively to study superconductivity in the 
Hubbard and t — J models^—. Go contains both normal 
(particle-hole) as well as anomalous (particle-particle) 
components in order to include superconducting pair- 
ing correlations. The Nambu spinor is defined by ^ d = 

(<£L, • ■ ■ , d\^, dii, • ■ ■ , G?4j,) and /i, v label the degrees of 
freedom within the cluster. 

Using a starting guess for the Weiss field Go, the clus- 
ter Green function G' is computed by solving a cluster 
impurity problem using a Lanczos exact diagonalization 
scheme, the details of which are discussed in Refs [23| or 
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with G' = - (Td^(r)4aM> and K» . = 
— (Td M -|-(r)d^4,(0)), the normal and anomalous time- 
ordered Green functions respectively. The cluster self- 
energy X' is obtained from 
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Finally, the following self consistency condition is em- 
ployed to recalculate Gq 1 iteratively until convergence is 



achieved, 

Go 1 far, 



dfc G(k, iui n ) 



E'(iw„),(4) 



with N C =A the cluster size and with the following defini- 
tion for the superlattice Green's function 



G(k,iu n ) 



fi - t(k) - S'(iw„) 
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where t(k) is the Fourier transform of the superlattice 
hopping matrix and the momentum integral is performed 
over the reduced Brillouin zone of the superlattice. 

The d-wave superconducting order parameter is de- 
fined as the expectation value of a particular pairing op- 
erator. Since the sites on the plaquette are connected 
via two distinct types of links (with hopping amplitudes 
t and t'), we define correspondingly two singlet pairing 
operators D and D' as 
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where the matrices D and D' are defined as follows: 
Dij = ±1 on t- finks in the x and y directions respectively, 
and likewise for DL on t' links. The corresponding order 
parameters are calculated using the anomalous part F of 
the superlattice Green function ([5]) as follows (for details, 
please consult Ref. [H): 



(7) 



and likewise for D' . The effective order parameter ^ is 
just the average (D + D')/2. 

The hopping strengths t and t' are defined keeping the 
average bandwidth to fixed to a constant value io of unity 

as 



t = tr 



At 



t' = t + At 
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with At > 0. At measures the degree of inhomogeneity 
in the system, which may be varied independently of the 
average bandwidth to- We insist on the importance of 
keeping the latter constant when varying At, since vary- 
ing to/U may cause effects that are likely more important 
that the inhomogeneity itself. 

An immediate question arises as to the appropriate 
choice for the cluster. As seen in Figure [I] there are 
three distinct ways in which a 2 x 2 plaquette might be 
selected, namely a plaquette with all t links; one with t 
links along x(y) axis and t' links along y(x) axis; or one 
with all t' links. We label them as plaquette A, B and 
C respectively. In order to select the appropriate clus- 
ter for the CDMFT calculation, one needs to determine 
which of the plaquettes best captures the physics of inho- 
mogeneity. Plaquette B seems to be the natural choice 
under such considerations, since it includes both t and 
t' links within the cluster allowing both to be treated 
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FIG. 2. (Color online) d— wave SC order parameters D,D' for the homogeneous case (t = t') for U = Sto, as a function of 
electron filling for plaquette A, B and C. The order parameter ^ = (D + D')/2 that is averaged over plaquettes A and C has 
been plotted alongside D,D' for plaquette B (middle plot) for comparison. In the latter case the curves overlap. 



on the same footing within CDMFT. In contrast, pla- 
quette A(C) treats only t(t') hoppings exactly while t'(t) 
are treated in a mean-field approximation via the bath 
degrees of freedom. 

Since the link inhomogeneity expands the unit-cell by 
a factor of two in each direction, the four-site plaquette 
really constitutes a single unit cell of the model, i.e., it is 
not really a cluster. As such, we are literally using single- 
site dynamical mean field theory instead of CDMFT, al- 
beit with a four-band model. Our treatment naturally 
collapses into a single band model (and CDMFT) in the 
homogeneous limit t! — t. 

In order to further justify our choice of cluster, we 
compute the superconducting order parameters {D,D') 
for the three clusters in the homogeneous limit (t = t'), 
where D and D' are expected to be identical. We see in 
Figure [5] that, whereas the order parameters D and D' 
are identical for plaquette B, they are significantly dif- 
ferent for plaquettes A and C. To understand this result, 
one has to remember that even in the homogeneous case, 
where t — t', the CDMFT lattice Green's function breaks 
translational symmetry, unless it is "periodized"— . As 
a consequence, the value of the order parameter on t 
and t' links will be different depending on whether they 
are inside or outside the plaquette. Mathematically, the 
self-energy matrix entering the lattice Green's function 
Eq.© is the same for all three plaquette choices when 
t = t'. However the D(k) operators in Eq.© have in- 
dices that are shifted with respect to those of the hopping 
matrix in Eq. ([5]) depending on which plaquette is chosen. 
This also explains why, for t = t', the value of D on pla- 
quette A is the same as the value of D' on plaquette 
C while for plaquette B, D = D' . It should be noted 
however, that, as illustrated by Fig. [5J the mean value 
of the order parameter averaged over the dissimilar links 
is virtually identical for any choice of plaquette in the 
homogeneous limit. 

Away from the homogeneous limit, we observe quan- 
titative differences for different choices of plaquette. In 
order to check whether or not the choice of cluster leads 
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FIG. 3. (Color online)Average d-wave superconducting or- 
der parameter as a function of electron filling at U = Sto, 
At = 0.05(top) and At = 0.1 0(bottom), for plaquette B, the 
average of results for plaquettes A and C, and the homoge- 
neous case. 



to a qualitative difference in our results, we compare the 
(i-wave order parameter for plaquette B to that of the 
average of plaquettes A and C for moderate levels of in- 
homogeneity (At = 0.05io and O.lio), as shown in Figure 
[3J We find that even though the average of the super- 
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FIG. 4. (Color online) (Top) The checkerboard Hubbard lad- 
der. (Bottom) d-wave order parameter vs t' jt for two different 
values of electron density n. Note that contrary to the rest of 
this paper, in this subsection we work directly with the ratio 
t' /t instead of At to have a fair comparison with published 
results. 

conducting order parameters for plaquettes A and C is 
larger than that for plaquette B, it never exceeds the cor- 
responding values for the homogeneous case. The results 
discussed in subsequent parts of this paper have been ob- 
tained using plaquette B as the cluster of choice. Note 
that since we do not compute direction dependent quan- 
tities, the results are identical for plaquette B, whichever 
of the two possible n/2 related orientations we choose. 



III. BENCHMARK WITH CHECKERBOARD 
HUBBARD LADDER 

To test the reliability of our CDMFT approach for in- 
homogeneous systems, we present results for the checker- 
board Hubbard ladder, in which the Hubbard model is 
denned on a one-dimensional, period two array of square 
plaquettes consisting of t links, connected by t' links as 
illustrated in the top panel of Figure [4] This is moti- 
vated by the availability of numerically exact results for 
this problem that were obtained using the Density Ma- 
trix Renormalization Group (DMRG) technique^. An 
enhancement in superconductivity, as measured by the 
pair binding energy, was found with increasing inhomo- 
geneity up to a moderately large value (t 1 /t ~ 0.6 for 
U=8t and tt,=0.875), where the pair binding energy at- 
tains its maximum value. While DMRG does not lead 
to long-range order on the ladder, the superconducting 
correlations decay algebraically. This will be mimicked 
by true long-range order in CDMFT, which treats long- 
range correlations in a mean-field way. 

Our studies on the checkerboard ladder find that the 
dependence of superconductivity on inhomogeneity is 
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FIG. 5. (Color online) Single-particle density of states in the 
superconducting state of the checkerboard Hubbard ladder 
for different t' jt at two values of the electron density n. 



qualitatively similar to the DMRG results. As seen in 
the bottom panel of Figure 21 the d-wave order pa- 
rameter increases with inhomogeneity (decreasing t' /t), 
and attains a maximum value at a relatively larger value 
of inhomogeneity {i! jt w 0.35). We do not have access 
to the pair binding energy, as opposed to DMRG, how- 
ever we can plot the density of states. The heights of 
the coherence peaks on each side of the energy gap in 
the one-particle density of states can be adopted as a 
measure of the strength of the superconducting correla- 
tions. These heights are indeed correlated with the mag- 
nitude of the order parameter on varying t' /t, as seen 
for two different dopings in Figure [5] Such behavior is 
in stark contrast to that of the two-dimensional checker- 
board Hubbard model, where, as we shall discuss below, 
there is no non-zero optimal value of inhomogeneity that 
favors superconductivity. Regardless of the contrasting 
observations in the two systems, this exercise serves to 
strengthen our claim of the validity of our CDMFT re- 
sults for the two-dimensional system. 
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FIG. 6. (Color online) d— wave SC order parameters D,D' 
and their average f as a function of electron filling n for U 
= 8to and At = 0.05 computed with plaquette B. 



IV. RESULTS 

This section is divided into three parts. We first de- 
scribe the results of computation of the d-wave supercon- 
ducting order parameter for strong and weak coupling, 
without allowing for antiferromagnetic long-range order. 
In addition, we show that in the superconducting state 
at a given doping, the heights of the peaks of the density 
of states lying on either side of the Fermi energy (across 
the energy gap) correlate with the magnitude of the or- 
der parameter. We then show that the correlation found 
previously^ between the low energy peak in the imagi- 
nary part of the spin susceptibility and the d-wave order 
parameter is still preserved in the inhomogeneous case. 
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FIG. 7. (Color online) (Top) d-wave superconducting order pa- 
rameter $ = (D + D')/2 as a function of chemical potential 
for various values of inhomogeneity At for plaquette B. (Bot- 
tom) Dependence of the lattice electron density on the chem- 
ical potential for the values of At that appear in the above 
plot. 



A. Superconducting order parameter 

Figure [6] shows D, D' and "F plotted as a function 
of the electron density n for U = 8to and At = 0.05. 
The strength of d— wave superconductivity over the entire 
doping range is larger across a link with larger hopping 
amplitude (t'). This is consistent with other studies 23 
which find that in the strong-coupling limit scales 
roughly with J = At 2 /U, the nearest-neighbor spin super- 
exchange coupling (at U = 8to we are indeed entering the 
strong-coupling regime). Quantitative differences aside, 
the plots of D and D' otherwise look qualitatively very 
similar, arising from the fact that t and t' are both treated 
on the same footing within the plaquette thereby allowing 
us to systematically isolate the physics of inhomogeneity 
from that of varying the effective bandwidth. 

The results for the superconducting order parameter 
for several values of inhomogeneity At, displayed in the 
top panel of Figure [3 exhibit two interesting features. 
We find that the strength of d-wave superconductivity 
decreases monotonically as a function of At over the en- 
tire doping range over which superconductivity exists. 



This stands in contrast to results obtained by studies 
on finite clusters^ - — where the pair binding energy was 
found to be maximized for moderately high levels of in- 
homogeneity at low hole dopings. However, our findings 
are in qualitative agreement with those of Doluweera et. 
aiiiS, in which DCA was used to study the checkerboard 
Hubbard model and the superconducting transition tem- 
perature T c was found to fall monotonically as a function 
of inhomogeneity in the underdoped regime. In addition, 
we observe a first-order superconducting to normal tran- 
sition in the underdoped regime beyond a moderately 
large level of inhomogeneity (At > O.I5to). The exis- 
tence of the first-order transition is confirmed by the ob- 
servation of hysteretic behavior in the order parameter 
depending on the initial state being normal (small dop- 
ing) or superconducting (larger doping) (Figure [8|), as 
the system is tuned across a superconducting to normal 
transition. Finally, as inhomogeneity is increased further 
to At > 0.16t (corresponding to t/t' < 0.72, we find that 
superconductivity is completely destroyed for all dopings. 
So far we have focused on strong coupling (U = 8to), 
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FIG. 8. (Color online) Average d-wave superconducting or- 
der parameter as a function of electron filling n for At= 0.15, 
showing hysteretic behavior as the system is tuned between 
the normal and the superconducting states in the underdoped 
regime. The blue (dashed) curve corresponds to increasing n, 
starting from a superconducting initial state close to optimal 
doping, while the red (bold) curve corresponds to reducing n, 
starting from a normal initial state close to half-filling. The 
value of the order parameter in the filling range of (0.89-0.94) 
is dependent on the initial state being normal or supercon- 
ducting. 



where the system is a Mott insulator at half- filling. It has 
been shown^ that in the strong-coupling regime, proxim- 
ity to the Mott insulating state leads to the suppression of 
superconducting order parameter close to half-filling. In 
contrast, the behavior is very different for weak coupling, 
where the Mott transition is absent, and no suppression 
is observed in the <i-wave order parameter in the under- 
doped regime unless antiferromagnetic long-range order 
is allowed^. Results for U — 4to in the presence of in- 
homogeneity are shown in Figure We find that the 
superconducting order parameter in the inhomogeneous 
case (At = O.lOio) is suppressed compared to the homo- 
geneous case except except for large dopings, where the 
superconducting order parameter is larger in the inho- 
mogeneous case. The maximum value of the order pa- 
rameter which occurs at half-filling in both cases, is how- 
ever larger in the homogeneous case. We must empha- 
size here that for weak coupling, the CDMFT method is 
not completely reliable as longer range antiferromagnetic 
correlations, which are important at weak coupling close 
to half-filling, are not adequately captured by this tech- 
nique. Nevertheless, these results may serve to demon- 
strate the qualitative difference between the weak cou- 
pling and strong coupling results. 



B. Density of states 

The principal disadvantage of using the d-wave or- 
der parameter as a measure of superconductivity is that 
the order parameter, though easily computed within 
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FIG. 9. (Color online) Average d-wave superconducting order 
parameter "5 = (D + D')/2, as a function of electron filling at 
U = 4tn, for homogeneous limit (At = 0) and At = O.lto for 
plaquette B. 



CDMFT, is not an experimentally measurable quantity, 
in contrast to the superconducting energy gap or T c . 
Within the BCS theory, the gap is given by the order 
parameter times the effective interaction^. Regardless 
of the validity of BCS theory for superconductivity in 
the Hubbard model, if one assumes that such a relation 
holds approximately true, a knowledge of the effective 
interaction would be required in addition to the order 
parameter to estimate the gap. In particular, in order 
to accurately estimate the gap as a function of At, one 
would have to determine the dependence of the effective 
interaction on inhomogeneity. 

A more straightforward and physically meaningful way 
to estimate the strength of superconducting correlations 
is to compute the single-particle density of states, which 
is directly measured in tunnelling experiments^, for 
varying At (Figure [TT)|) . The density of states fea- 
tures a gap around the Fermi energy as seen in earlier 
studies^.. However, the magnitude of the energy gap 
does not change appreciably with At on the scale of the 
Lorentzian broadening rj— O.lto, used to compute the 
density of states, which makes it unsuitable to probe re- 
liably the variation of superconductivity with inhomo- 
geneity. Alternatively, one might consider the quasipar- 
ticle peak heights in the density of states. The quasi- 
particle spectrum of BCS superconductors is character- 
ized by peaks on either side of the energy gap, and the 
height of the peaks is a measure of the coherence of the 
quasiparticle excitations, and may therefore be used to 
gauge the strength of coherence in the superconducting 
state for varying At at a fixed doping. As seen in Fig- 
ure [TUl for various At the heights of the peaks vary 
concomitantly with the magnitude of the d-wave order 
parameter (Figure [7]) in both underdoped and overdoped 
regimes, thereby providing an experimentally measurable 
probe of superconductivity whose behavior is consistent 
with that of the order parameter. 
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FIG. 10. (Color online) Single-particle density of states for 
different values of At at two values of the electron density n 
for plaquette B. 



FIG. 11. (Color online) Imaginary part of the cluster spin 
susceptibility at Q=(n,n) (top) and Q=(0, n) (bottom), for 
three different At values, with U =8to and n=0.84. 



C. Spin susceptibility 

The top panel in Figure [11] displays the imaginary part 
of the cluster spin susceptibility x"( w ) m the overdoped 
regime (n= 0.84), for U= 8to at Q = (tt, tt) for three dif- 
ferent values of At. The strength of the low energy peak, 
whose connection with superconductivity was confirmed 
by theoretical studie s 26 i 28 on the overdoped homogeneous 
Hubbard model and neutron scattering experiments 2 ^— 
on LSCO samples, falls with increasing inhomogeneity, 
concomitant with the behavior of the superconducting 
order parameter. In contrast with the homogeneous case, 
the fall in the susceptibility peak with inhomogeneity is 
slower than that of d-wave superconductivity. Neverthe- 
less, we find that the association between superconduc- 
tivity and the low energy peak in the spin susceptibility 
at Q ~(tt, tt) in the overdoped regime is preserved in the 
presence of inhomogeneity as well. As seen in Figure 
the low energy peak in the spin susceptibility is present 
in the superconducting state, while it disappears in the 
normal state ( 1 F = 0), which further corroborates the 
connection between antiferromagnetic fluctuations and 
d-wave superconductivity found previously in the homo- 
geneous Hubbard model. 



Another interesting feature in the spin susceptibility is 
the enhancement of the Q =(0,7r) component of x"( w ) 
with inhomogeneity as seen in the bottom panel in Fig- 
ure 111! which may be an indication of development of 
spin fluctuations competing with the predominant anti- 
ferromagnetic fluctuations in the system. One must how- 
ever be careful in interpreting the above results, since the 
cluster which is inherently anisotropic, favors such spin 
fluctuations, and does not necessarily reflect on the na- 
ture of the long-range spin correlations. The strength 
of the low energy peak (around w=0.3) however remains 
significantly smaller than that at Q = (7r,7r), indicating 
that antiferromagnetic correlations, although weakened 
by inhomogeneity, continue to dominate the physics of 
the two-dimensional Hubbard model in the presence of 
moderate checkerboard-type inhomogeneity. 

V. DISCUSSION 

Our studies indicate that in the strong coupling 
regime, which is considered relevant for the cuprates, 
checkerboard-type inhomogeneity on the square lattice 
Hubbard model is detrimental to d-wave superconduc- 
tivity over the entire doping range of interest. Supercon- 
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FIG. 12. (Color online) Imaginary part of the cluster spin 
susceptibility at Q—(n,ir) in the superconducting (At=0.12, 
n=0.84), and normal state (A£=0.12, n=0.76 and Af=0.20, 
n=0.77) for U=8t . 

ductivity is completely destroyed beyond a moderately 
large inhomogeneity level. This is a strikingly different 
conclusion from what was reported in some recent works 
that studied the highly underdoped regime of the checker- 
board Hubbard model on finite clusters using exact 
diagonalization 10-12 , and contractor renormalization 13 
methods, where the pair binding energy was found to 
be maximum at moderate levels of inhomogeneity at low 
hole concentrations. 

The difference between these results and those of the 
present study may be attributed to several reasons. 
Firstly, CDMFT is substantially better equipped in cap- 
turing the physics of the extended lattice, compared to fi- 
nite clusters which are expected to have significant finite- 
size effects. Therefore, interpreting the results obtained 
on finite systems and using them to predict the nature 
of superconductivity in the extended system must be un- 
dertaken cautiously. At least CDMFT captures some of 
the physics of the infinite lattice in a mean-field way, even 
though correlations are taken into account exactly only 
on short length scales. Secondly, previous studies com- 
puted the pair binding energy to quantify superconduc- 
tivity which, in contrast to the <i-wave order parameter 
considered in our study, is not a measure of supercon- 
ducting phase coherence in the system. Furthermore, 
our results are in qualitative agreement with those of 
Doluweera et al. 1 ® , who used DC A on a four-site pla- 
quette to study superconductivity on the checkerboard 
Hubbard model at finite temperature. Their study re- 
ported a monotonic suppression of T c as a function of 
inhomogeneity. In contrast, our work focused on the d- 
wave order parameter at T = as a measure of super- 
conducting strength, providing an alternative approach 
to this problem. It is noteworthy that in contrast to the 
aforementioned work where a single plaqucttc configura- 
tion was studied with uniform hopping within the cluster 
(corresponding to plaquette A(C) in our work depend- 
ing on t' being larger (smaller) that i), we have verified 



the dependence of results on the choice of cluster. It is 
not surprising, however, that our results are qualitatively 
similar, since DCA and CDMFT are both self-consistent 
cluster methods that effectively capture the short-range 
correlations of the system, and whose results become ex- 
act in the limit of infinite cluster size. 

There is a qualitative difference between the strong and 
weak coupling results in the extreme overdoped regime 
where there is an enhancement in superconductivity at 
weak coupling that is not observed at strong coupling. It 
is also worth mentioning here that our investigations in- 
clude the overdoped regime of the superconducting phase 
of the checkerboard Hubbard model, which has not been 
considered in any of the previously mentioned studies. 

In the strong coupling case, the gradual suppression of 
the order parameter with inhomogeneity, followed by a 
first-order transition to the normal state at moderately 
large inhomogeneity in the underdoped regime has not 
been noticed before. This result should be interpreted 
in light of a recent study that demonstrates that a first- 
order metal-metal transition lies beneath the supercon- 
ducting dome 3 - 2 - and that this transition is directly linked 
to the Mott transition. The superconducting phase of 
the unusual metal found close to half-filling is more sen- 
sitive to inhomogeneity. This can be verified by looking 
at Fig. [7] at fixed filling as a function of inhomogeneity: 
the figure clearly suggests that a first order transition to 
the normal state occurs. 

Finally, our results support the connection between an- 
tiferromagnetic fluctuations and superconductivity found 
previously for the homogeneous case i 26 ' 33 Indeed, the cor- 
relation between the magnitude of the superconducting 
order parameter and the height of the first peak in \" in 
the overdoped regime 2 ^ remains valid for the inhomoge- 
neous case studied here, although the fall of the suscep- 
tibility peak is slower than that of the order parameter. 
It is entirely plausible therefore that a suppression in the 
superconductivity is tied to the of weakening of antiferro- 
magnetic correlations in the presence of a checkerboard- 
type inhomogeneity in nearest-neighbor hopping. 



VI. CONCLUSION 

We benchmarked our approach with numerically exact 
DMRG results on the checkerboard Hubbard ladder 25 . 
The quantities that can be obtained in CDMFT and in 
DMRG are different but qualitatively both approaches 
show similar results, namely that there exists an optimal 
inhomogeneity for superconductivity on a ladder. 

Despite suggestions that this result is general, namely 
that there is always an optimal inhomogeneity for 
superconductivity ^ 15 ' 34 " —, our study demonstrates in- 
stead that this statement may not be valid for d-wave su- 
perconductivity in the two-dimensional Hubbard model 
in the presence of arbitrary types of inhomogeneity. In- 
deed, previous DCA results at finite temperature^ and 
our CDMFT study at zero temperature both find that 
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when inhomogeneity on the checkerboard Hubbard model 
is in the nearest-neighbor hopping, either the maximum 
T<~ or the maximum value of the superconducting or- 
der parameter at zero temperature cannot exceed that 
of the homogeneous system. The size of the order pa- 
rameter at a given doping can be taken as a measure of 
the strength of superconductivity, since we have shown 
that it manifests itself directly in the height of the co- 
herence peaks in the density of states. Note that since 
one finds, with CDMFT, that site inhomogeneity on the 
checkerboard lattice does lead to an optimal inhomogene- 
ity for superconductivity,- 15 quantum cluster methods do 
not have intrinsic limitations that prohibit finding en- 
hanced superconductivity in the presence of inhomogene- 
ity. 

For the model of interest, we have explored a larger 
doping range than previous studies as well as both weak 
and strong coupling regimes. In the weak coupling case, 
as a result of inhomogeneity, superconductivity is sup- 
pressed in the underdoped and enhanced in the over- 
doped regime, which however does not surpass the max- 
imum possible value of the order parameter in the ho- 
mogeneous case. In the strong coupling case, our results 
can be summarized by the following observations a) a 
monotonic suppression of the superconducting order pa- 
rameter for all dopings and b) at given inhomogeneity, a 



first order transition between normal and superconduct- 
ing state at finite doping on the underdoped side, with 
the superconducting dome disappearing suddenly on fur- 
ther increasing the inhomogeneity. 

Further research could look into the effect of next- 
nearest-neighbor hopping or of mixed types of inhomo- 
geneity, including both hopping and site energies. This 
would help again to verify the generality of the connec- 
tion confirmed here between antiferromagnetic fluctua- 
tions and superconductivity . 26 i 33 One should also explore 
more closely the relationship between the domain where 
a first order transition is induced by inhomogeneity at 
strong coupling in the underdoped regime and the do- 
main where a first order transition between two metals 
was found recently for the homogeneous casej 3 ^ The lat- 
ter phenomenon was clearly linked to the Mott transition. 
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